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Abstract 

We define a family of asymmetric processes for particles on a one-dimensional lat- 
tice, depending on a continuous parameter A 6 [0,1], interpolating between the com- 
pletely asymmetric processes [1] (for A = 1 ) and the n = 1 drop-push models [2] 
(for A = ). For arbitrary A, the model describes an exclusion process, in which a 
particle pushes its right neighbouring particles to the right, with rates depending on 
the number of these particles. Using the Bethe ansatz, we obtain the exact solution of 
the master equation . 
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1 Introduction 



Various versions of one-dimensional asymmetric simple exclusion processes (ASEP) have 
been shown to be of physical interest in a variety of problems including the Kinetics of 
biopolymerization [3], polymers in random media, dynamical models of interface growth [4], 
and traffic models [5]. This model is also related to the noisy Burgers equation [6], and thus 
to the study of shocks [7,8]. Besides the equilibrium properties of this model, its dynamical 
properties have also been studied in [8-10]. 

Recently the totally ASEP model, with sequencial updating on an infinite lattice, has 
been solved exactly by G. M. Schiitz [1] using the coordinate Bethe ansatz. In this model, 
each lattice site can be occupied by at most one particle and a particle hops with rate one 
to its right neighbouring site if it is not already occupied; otherwise the attempted move is 
rejected. In his work, instead of using the quantum Hamiltonian formalism, which is suit- 
able for studying the dynamical exponents and certain time-dependent correlation functions, 
Schiitz adopted the coordinate representation for writing the master equation. By solving 
the master equation exactly, he was able to obtain explicit expressions for conditional prob- 
abilities P(xi, #2, • • • , %N] t\yi, 2/2, • • • , UN] 0) of finding N particles on lattice sites xi, • ■ ■ xn 
at time t with initial occupation yi, • • • , yjv at time t — 0. 

The master equation for the probability of finding particle 1 on site k\ , particle 2 on site 
k 2 , ■ • • , and particle N on site k N , with k N > k^-i > ...k 2 > hi, is written as 

d 

— P(ki, k 2 , ■ • • , k N , t) = P(k x - 1, k 2 , • ■ • , k N , t) + P(h, k 2 - 1, • • • , k N , t) H 

+ P(k u k 2 , ■ ■ ■ , k N - 1, t) - NP(k u fc 2 , • • • , k N , t), (1) 
if ki + \ — ki > 1. This equation was then augmented by the following boundary condition 

P(k,k,t) = P(k,k + l,t), VJfe. (2) 

In writing (2), we have supressed for simplicity the position of all the other particles, bearing 
in mind that this condition should hold for every pair of adjacent variables ki and k i+ i. In 
the following we always use this simplified notation. It was then shown that (1) and (2) give 
the correct master equation in the whole physical region ( i.e. the region ki < k i+ i ) for the 
probabilities. In the rest of [1], the exact solution of the master equation (1) (with boundary 
condition (2)) is constructed. 

We now describe what we have done in the present paper. In section 2, we substitute 
the boundary condition (2) by 

P(k,k,t) = P(k-l,k,t), Vk, (3) 

and show that this boundary condition, together with (1), describes the n — 1 drop-push 
dynamics [2]. In this process, even if the right neighbouring sites of a particle are occupied, 
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the particle hops with rate one to the next right site, pushing the right neighbouring particles 
to the next sites. This means that all the following processes occur with equal rate one : 



AO -> OA, 

AAO -> OAA, 
AAAO -> OAAA, 

AA---A0 -> 0AA---A (4) 

n n 

where we have adopted the standard notation for representing a particle by A and a vacancy 
by 0. We then obtain a closed form for the conditional proabilities for this process. 

This process, in which a particle pushes as many particles with rate one, is the opposite 
extreme of what was solved by Schiitz, and interestingly admits a closed form solution for 
the conditional probabilities P(xi,x 2 , ■ ■ ■ , x^; t\yi, y 2 , - ■ ■ , Vn] 0) in the form of an N x N 
determinant. 

In section 3, we combine the boundary condition (2) and (3) in the form 

P(k,k,t) = \P(k,k + l,t) + (l-\)P(k-l,k,t), Vfc, (5) 

and show that the resulting master equation ((1) and (5) ), describes a processin which the 
processes shown in (4) occur with unequal rates: namely the process 

AAA---A0 -> 0AAA---A, (6) 

n n 

occurs with rate 

(7) 



1 + A + ( A )2 + ... + ( A)n> 



where fj, — 1 — A. We call this model Generalized Totally Asymmetric Exclusion Process. In 
the limit A — > 0, we have r n = 1, Vn, and in the limit A = 1, we have r = 1 and r n ^ = 0. 
Note also that r n+ i < r n , Wn. Therefore this process is perhaps more physical than the two 
extreme cases studied in [1] and in section 2 of this paper. 

In section 4, we use the coordinate Bethe ansatz and solve the master equation of the 
process defined in section 3, and show that there is no bound state in the spectrum. 

In section 5, we write the quantum Hamiltonian formalism for the generalized process 
and, using a particle-hole exchange transformation, show that this generalized process is 
equivalent (i.e. in the same universality class ) to another process, where particles hop only 
to the left. In this new process, if a left neighbouring site is occupied, the move is rejected, 
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but if a set of the left neighbouring sites are empty, the particle hops with distance dependent 
rates to these sites: 

OA -> AO with rate r 

00A -> A00 with rate r x 

000---0A^ A000---0 with rate r„. (8) 

ra n 

Therefore a transformation as simple as a particle-hole exchange, when applied to our gen- 
eralized process, has an interesting physical consequence. Models with different values of 
A, all allow exact solutions in the form of coordinate Bethe ansatz and their spectrum have 
only the continuous part, but only the limiting cases of these models (A = and 1 ) allow a 
closed solution in the form of a determinant. 
We end up the paper with conclusion in section 5. 

2 Generalized totally asymmetric exclusion process 
with A = 

We augment the master equation (1) with the boundary condition (3). Although we derive 
the rates for arbitrary A in the next section by a general argument, here we want to show 
that for A = case, the master equation (1) (together with the boundary condition (3)) 
describe an n = 1 drop-push dynamics. For simplicity, consider the two particle sector of 
n = 1 drop push-dynamics. The master equations are 

^P(k 1 M,t)=P(k 1 -l,k 2 ,t) + P(k 1 ,k2-l,t)-2P(k 1 ,k 2 ,t), k 2 >h + l, (9) 
d 

—P{k, k + l,t)=P(k-l,k + 1, t) + P(k - 1, k, t) - 2P(k, k + l,t). (10) 
at 

Now, if we choose the boundary condition 

P(k,k,t) = P(k-l,k,t), VJfe, (3) 

eq. (10) can be written as 
d 

—P{k, k + l,t) = P(k-l,k + 1, t) + P{k, k, t) - 2P(k, k + l,t), (11) 
at 

which is of the same form as (9). 

In the three particle sector, the extra equation which needs to be taken into account is 
d 

— P{k, k+1, k+2) = P(k-1, k+1, k+2)+P(k-l, k, k+2)+P(k-l, k, k+l)-3P(k, k+1, k+2). 

(12) 
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Using the boundary condition (3), the second and the third terms on the r.h.s. of (12) can 
be written as 

P(k-l,k,k + 2) = P(k,k,k + 2), (13) 
P(k - 1, k, k + 1) = P{k, k, k + 1) = P(k, k + 1, k + 1), (14) 
which means that eq. (12) is equivalent to the following standard form 
d 



dt 



P(k, k+1, k+2) = P(k-1, k+1, A;+2)+P(A;, k, k+2)+P(k, k+1, k+l)-W{k, k+1, k+2). 

(15) 

This procedure can be repeated for any sector. We will give a general proof in the next 
section. 

To solve the master equation (1), and the boundary condition (3), for the conditional 
probability P(xi, x 2 , ■ • • , xjy; t\yi, y 2 ,---, Vn] 0), we set, following Schiitz [1], 



P(x u x 2 , ■ ■ ■ , x N ; y 2 ,---, Vn] 0) = e Nt detG N , 



(16) 



where is an iV x matrix with entries G^j — gi—jixi — yj, t). The functions g p (x,t) are 
to be determined. Writing Gn as 



where 



" G 1 (x 1 ,t) 
G N = det Gi(xi, t) 
_ G N (x N ,t) 

Gi(xi) = (gi-^Xi -yi,t),gi- 2 (xi - y 2 ,t),- • • ,gi- N {xi -y N ,t)), 



(17) 



(18) 



N 

Edet 



i=i 



obtain 






" G x (x u t) 






-Q^Gi^Xi, t) 


N 

= £ det 

i=i 


Gi(xi - l,t) 


Gn{xni t) 




GN{XN,t) 



(19) 



the solution of which is 



d_ 

dt 



Gi(xi,t) = Gi(xi - l,t). 



(20) 
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in the boundary condition (3), we 


obtain 




"l\, J 'l) 1 ) 




Gi(x!,t) 


det 


G k -i{x,t) 


= det 


G k -i(x - l,t) 




G k (x,t) 




G k (x,t) 




_ G N (x N ,t) _ 




G N (x N ,t) 



(21) 



the solution of which is 



G k -!(x, t) = G k -x{x -l,t) + pG k (x, t), 



(22) 



where (3 is an arbitrary parameter. The explicit form of the function g p (x,t) can now be 
determined: these functions, as seen by eqs. (20) and (22), should satisfy the following 
relations 



d 



dt 



9 P (n,t) = g p (n - l,t), 



(23) 

g p (n,t) = g p (n- l,t) + (3g P+ i(n, t). (24) 

Defining the generating functions (or ^-transforms) g p (z,t) := X^_oo ^Qvi 71 -,^)-, e Q s - (23) 
and (24) are converted to 

d _ 

—g p (z,t) = zg p (z,t), (25) 

and 

9p+i(z,t) = -(1 - z)g p (z,t), (26) 
the solution of which is simply obtained as 

.1-2, 



g p (z,t) = e z %(z,0) = e zt (—) p g o (z,0). 



(27) 



go(z, 0) is nothing but the generating function for go{n, 0), the one particle sector probabilities 
at t — 0. Since P(x, 0\y, 0) = go(x—y, 0) = 5 XtV , we have go(n, 0) = 5 nt o, and thus go(z, 0) = 1, 
giving finally 

I - • 

(28) 



~g P ^t)=e z \ l —^y. 



The parameter f3, as long as it is nonzero, drops out of the determinant and so we can set it 
equal to unity. The functions g p (n, t) are obtained by expanding the generating functions. 
Note that the functions g p (z, t) should be expanded in terms of positive powers of z, if p < 0. 
This is due to the fact that, for p < 0, as n — > — oo, the function g p (n,t) tend to zero, since 
this limit is in the physical region. This expansion yields, formally, 



(29) 
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If p > 0, g p (n, t) is converted to a finite sum 

kj (n — k)\ 



9 P >o(n,t)= £ PUJ t »-*. (30) 



fc=0 

In particular, 



<7o(M) = ^- (31) 



If p < 0, g p (n, t) is converted to another finite sum 

* (TM)= £l ipi-i (32) 

We have thus obtained an explicit relation for the conditional probability. 



3 Generalized totally asymmetric exclusion process 
with arbitrary A 

We now consider the master equation (1) together with the boundary condition 

P(k,k,t)=\P(k,k + l,t)+fjiP(k-l,k,t), Vife. (5) 

It can be easily shown that the conservation of probability demands that fj, — 1 — A. In order 
to understand what type of process is described by these equations, we first look at the two 
particle case. Eqs. (1) and (5) yield 

^P(k,k + 1) = P(k-l,k + l) + P(k,k)-2P(k,k + l) 

= P(k-l,k + l)+ ^P(k-l,k) - (1 + fx)P(k,k + l), (33) 

which means the following rates. 

AO -> OA with rate r = 1, 

AAO -> OAA with rate n = /i. 

To find the rates in the general case, we first prove a lemma. 
Lemma: Equation (5) implies, for arbitrary n, the following 

P(k, k+1, k+2, • • • , k+n-1, k+n, k+n) = (l-r n+1 )P(k, fc+1, A;+2, • • • , k+n-1, k+n, k+n+1) 
+ r n+1 P(k -l,k,k + l, - ■ ■ ,k + n- 2,k + n-l,k + n), (34) 

where 

rn = (l + ^ + A 2 + -- + Ar 1 - (35) 
6 



Proof: We proceed by induction. For n — 0, eqs. (34) and (35) are the same as (5), as 
r 1 = ji. Assuming now, that (34) and (35) are correct for n — 1, and using eq. (5), we have 

P(k, k + 1, • • • , k + n — 1, k + n, k + n) 

= XP(k, k + 1, - ■ ■ ,k + n-l,k + n,k + n + l) + nP{k, k + 1, • • • ,k + n - l,k + n - l,k + n) 

= XP(k, k + l,---,k + n-l,k + n,k + n + l) 

+/i{(l — r n )P(k, k + 1, • • • , k+n — 1, k+n, k+n)+r n P(k — l, k, • • • ,k+n — 2,k+n — l,k+n)}, 

(36) 

or 

P(k, k + 1, • • • , k + n — 1, k + n, k + n) = s n+i P(k, k + 1, • • • ,k + n — l,k + n, k + n + 1) 
+ r n+1 P(k — l,k,---,k + n — 2,k + n — l,k + n), (37) 

where 

A fir n _ , . 

jz r — s n +i , z 7- r — r n+1 . [66 

l-fjL{l-r n ) l-n(l-r n ) 

From eq.(38), it is seen that s n+ i +r n+ i = 1. One can now solve the second equation of (38) 

for r n+ i to obtain 

which gives 



^ 1 = i + -{i + - + (-) 2 + --- + (-n = i + - + (-) 2 + --- + (-r 1 - 



/i \X /J, pL 11 fi fi 

This proves the lemma. 

We now consider a collection of n adjacent particles and write the master equation for 
this configuration by eq. (1): 

dP 

—(k, k + l,k + 2,---,k + n-l) = 
at 



n— 1 



P(K k+1, • • • , k+i-2, k+i-1, k+i-1, k+i+1, • • • , k+n-l)-nP(k, k+1, k+2, • • • , k+n-l). 

i=0 

(39) 

Using (34), we find 

dP n ~ x 

— (k, k + 1, k + 2, k+n-l) = V nP(k-l, k,---, k+i-2, k+i-1, k+i + 1, ■■■ , k+n-l) 

~ (E r ^ k + 1, k + 2, k + n-l). (40) 

It is now obvious that the above equation describes a process in which a collection of % + 1 
adjacent particles hop to the right with rate r\, as claimed in the introdution. 
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4 The Bethe ansatz solution for arbitary A 

In this section we denote the position of the particles by Xi G Z rather than ki, and apply 
the Bethe ansatz for the solution of the master equation (1) and the boundary condition (5). 
Writing P N (x±, • • • , xn, t) = e~ ejv *\l/ n(%i, • • • , xn), will turn (1) into an eigenvalue equation 
for ^ N (x!, ■ ■ -,x N ): 

N 

E Vn(xi, •'•,x i -l,---,x N ) = (N- e N )^ N (x 1 , ■ ■ ■ , x h ■ ■ ■ , x N ). (41) 

We write the coordinate Bethe ansatz for \1> in the form: 

vM*i,"--,^) = E^e^- x , (42) 

where x and p stand for the n-tuple coordinates and momenta and cr(p) is a permutation 
of momenta. The sum is over all permutations. Inserting (42) into (41) yields 

E A a e ia ^- X (e' la ^ + e- taip ^ + ■■■ + e" iCT(pjv) ) = (N - e N )^ N {x 1 , ■ ■ ■ , x N ). (43) 

c 

The sum in the paranthesis can be written as J2k=i e~ Wk and taken outside yielding 

N N 

e J v:=E^) = E( 1 -^ Pfc ). (44) 



k=l k=l 



Note that due to translational invariance, is also an eigenvector of total momentum P, 
which in the lattice is defined as the logarithm of the shift operator U = e~ lP : 

(U^ N )(x ir --,x N ) :=V N (xi-l,x 2 -l,---,x N -l). (45) 

Acting by U on (42), we obtain 

(PV N )(x 1 ,---,x N ) = (p 1 + ---+p N )^ N (x 1 ,---,x N ). (46) 

Therefore the eigenvectors have additive total energies and momenta. Inserting (42) in 
the boundary condition (5), rewritten in an unabbreviated form: 

(xi, • • • , Xi = f , x i+1 = £, • • • , x N ) = A*(xi, • • • , Xi = f , = £ + 1, • • • , %n) 

+H i &(x 1 , ■ ■ • , Xi = £ - 1, a; i+1 = f , • • • , x N ), 

we obtain 

^^E^i.i+i^fe^+'^^O+^CPi+i)^ ^ _ ^ e ia( Pl+1 ) _ ^g-^CPijy = o. (47) 
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We denote the expression in the bracket by B c . Noting that the prefactor is unaffected by 
an interchange of Pi and Pi+i, it follows that the proper coefficient of each prefactor, which 
should vanish, is B a + B aa . , where <7j is the generator of Sn (the permutation group of N 
object ) which only interchanges pi and pi + i 

a i (p 1 ,---,p i ,p i+l ,---,p N ) = (pi,---,Pi+i,Pi,---,p N ), (48) 

and cCj stands for the product of two group elements, a acting after <jj. Therefore we find: 

4,(1 - \e m{ P^ - fie- ia{p ^) + A ac .{\ - \e ia ^ - /le"^^) = 0, 

or 

= + fie _i = <j(p i+1 )). (49) 

This relation, in effect, allows one to find all the A^s in terms of Ai (which is set to unity). 
The first few coefficients, corresponding to the elements 1, a±, a 2 , <7i<72, o<iO\, <j\<7i<7\ are : 

A\ — 1 , A Cl = Su , A a2 = 5*23, 

A aia2 = S12S13 , A U2<Tl — S13S23 , A aia2(Tl = S12S13S23, (50) 

where Sij = S(pi,pj). The form of the scattering matrix Sy could also be found from the 
two particle sector alone. The above analysis shows in fact the factorizibility of the S matrix 
in the general case, a sign of the integrabilty of the problem. 
To find the range of p^s, we analyze the S matrix 

Ae^+/ie-^-l C2i 

^ 12 = ; r~ ~ — — > (51) 

1 - Ae^ 1 - yue"^ 2 C12 

and the two particle wave function 

* 2 ( Xl ,x 2 ) = Cl2 e i(piiE1+P2a:2) + C2ie <(p2iei+pia!2) , 

or 

x) = e lPX (c 12 e ipx + c^e"**) , (52) 

where X := \{x\ + x 2 ), x := X\ — x 2 , P '■— P\ + £>2 and p := \{p\ — P2) with clear physical 
meanings. Since x is negative ( x\ < X2 ), to have a bound state one of the following set of 
condintions should be satisfied simultaneously. Either 

C12 = , Imp > , ImP = 0, 

or 

c 2 i = , Imp < , ImP = 0. 
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Rewriting c 12 m terms of the new momenta, we find 

c 12 = e ip (e- tp - \e 1 ^ - fie~^). 

Since ImP = we have 

|Ae*T + ne~ i % \ < X + /jl = 1. 



(53) 



(54) 



Noting that |e~* p | > 1, it is seen that C12 can not vanish. A similar analysis applies for the 
second set of conditions. Therefore no bound state exists in the spectrum, and the range of 
all momentum variables is [0, 27r). 

To find the conditional probability P^(xi, x 2 , x^', t\yi, y 2 , ■■■Un] 0), one should take a 
linear combination of the eigenf unctions ^n, with suitable coefficients. Consider the two 
particle sector. We have 

P 2 (x 1 ,x 2 ;t\y 1 ,y 2 ;0) = F F ^^ e -^^\^^^^ x ^ X2 y (55) 

JO JO 2,71 LlX 

This is just a linear combination of the eigenfunctions, satisfying the initial condition 



P 2 (x!,x 2 ] 0|yi, y 2 ; 0) = S xuyi 5 X2; y 2 , 

in the physical region (x 2 > x±, y 2 > yi). The eigenfunction ty 2 (xi,x 2 ) in (55) is normalized 
according to 

V 2 ( Xl ,x 2 ) = e iipiXl+P2X2) + S 12 e l(p2Xl+PlX2) . 

To avoid the singularity in Si 2 , we set p± — > pi + ie. With this prescription, the contribution 
of the second term in \I/ 2 to P 2 (xi, x 2 , 0\yi, y 2 , 0) identically vanishes in the physical region. 
Using the variables £ := e rpi and rj := e~' ip2 , a simple contour integration yields 



P 2 (x 1 ,x 2 ;t\y 1 ,y 2 ;0) 



-2t 



-ffi-yi 



(xi -yi) \ (x 2 -y 2 )\ 



oo k 



j-X\—y2—k+m 



U\ +X2-yi+rn 

_ I I \ m ^, k ~ m 

\m y (x 2 — yi + m)! (xi — y 2 — k + m)! 



fc=0 m=0 



X 



At 



^1-^/2-^ + ^+1 



p{x2 -yi + m) 
t 



(56) 



It is easy to see, explicitly, that this solution satisfies the initial condition in the physical 
region. Also, in the limiting cases A = 1 and A = 0, it reduces, respectively, to 



P2(x 1 ,x 2 ;t\y 1 ,y 2 ;0) = e 



-2t 



fxi-yi 



fX2-y2 



. (xi - yi)\ (x2 - 1/2)! 



[x 1 -y 2 )\ (xi - 2/2 + I)! 



00 j-x2-yi+k \ 

Sfe -„ + *).}' (57) 



10 



obtained in [1], and 

P 2 (xi,x 2 ;t\yi,y 2 ,0) = e~ 2t 



. (xi - yi)\ (x 2 - V2V- 

4.X2—U1 fX2—yi — l 



(x 2 -yi)\ (x 2 -y 1 -l)\ 

obtained in the present paper. 

The treatment of the N particle case is similar. We have 



, (58) 



P N ( Xl , -~,x N ;t\ yi ,-~, y N ; 0) = T ^ • • • (** ^ e -E«'-EM. %(ll) . . . , Xn) . 

Jo Ztc Jo Z7T 



The integration is defined with the following e-prescription: in S^- (i < j), pi is replaced by 
Pi + it- 



5 Hamiltonian approach 

The Hilbert space of generalized totally asymmetric exclusion process is H = ®C 2 , the 
tensor product of all the local Hilbert spaces of the lattice sites. C 2 is the two dimensional 

vector space with basis states >= [ J j and |1 >= ( ^ I '-1 |(> states K' > an( l > 



represent vaccant and occupied sites, respectively. The local operators rii 




1 



<jT = \ I , and (T- — [ I are the number , annihilation, and creation operators, 

1 o o ) V 1 

respectively. Their action on a bra state < a|, (a = 0, 1), can be conveniently represented 
as < a\n = a < a\, < a\a + — (1 — a) < 1 — a\ and < a\a~ = a < 1 — a\ . 

The state of the system \^f(t) > evolves according to the Schrodinger type equation 
— ■^\^(t) >= H\^(t) >. The connection between the two representations is given by the 
relation 

P(h,k 2 ,---,k N ,t) =< k 1 ,k 2 ,---,h N \V N (t) >=< 0\a+a+-- -,a+ N \^ N (t) > . (59) 
The Hamiltonian of the process can now be written as 

ff = -EE r i-iW0-^))> (eo) 

fceL i>i 

where L represents the sites of the lattice and 

v k{l) = (?t n k+ink+2 ■ ■ ■ n k+ i-ia^ +h (61) 
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w k (l) = n k n k+l n k+2 ■ ■ -n fc+ i_i(l - n k+l ). 



(62) 



Consider a bra state containing / particles on adjacent sites: < k + 1, k + 2, • • • , k + l\. The 
only terms in H with nonvanishing action on this state are 

< k + l,k + 2,---,k + l\v k (i) =< k,k + l, - • • ,k + i — l,k + i + l, - • • ,k + l\, l<i<l, (63) 

< k + l,k + 2, ■■■,k + l\w k+1+l _ i (i) =< k + l,k + 2, •••,fc + Z|, l<i<Z. (64) 

Note that the action of the above operators on every other state which contains, beside the 
above particles, other collection of particles disconnected from the above one is the same. 
Using eqs. (59) to (64), one arrives at eq. (40) for the evolution of the probability. Note that 
the quantum Hamiltonian (60) is a stochastic operator, meaning that all of it's off-diagonal 
matrix elements are non-positive with the sum of entries in each column being equal to zero. 
This last property is expressed by saying that < S\H — where < S\ is the sum of all basis 
states of Ti. Equivalent models may be obtained by constructing operators Q : 7i — > 7i and 
Hamitonians H' = QHQ^ 1 , which preserve the above properties. An obvious example is the 
particle-hole exchange operator = Ili a f- It clearly has the property that < S\Q =< S\, so 
that for H' we also have < S\H' — 0. It is easy to see that the this transformation induces 
the changes n <-> 1 — n and <j + <-> o~ . So the master equation obtained from H' describes 
the process (8). 



We have defined a generalized exclusion process, parameterized by a real parameter A G [0, 1], 
and have shown that the master equation of this model admits for every A an exact solution 
via the coordinate Bethe ansatz. We have also shown that this model interpolates continously 
between two very different models: the totally asymmetric exclusion model (for A = 1), which 
we may consider as the weakcoupling limit and the drop-push model (for A = 0), which may 
be considered as the strong coupling limit of the model. In these two limits, the solution 
acquires a simple determinant form. 

Our work can be further investigated in one definite way. It may be that the point A = \ 
is a point of phase transition and the study of the equilibrium properties of the model on 
a periodic lattice may reaveal this transition. There are already two pieces of evidence for 
the validity of this conjecture. First, there is some sort of duality between two models two 
models with parameters symmetric with respect to A = 1/2. To be more specific, we have 
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S(pi,P2] A) = S(-p 2 , -pi] 1 - A). 



(65) 
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Second, the large I behaviour of the transition rates is 




1 X — 1 

T> A - 2 



It will be interesting to study the stationary behaviour of this system along the lines which 
have been followed in [11-14], to see what kind of phases develop in the system by varying 
A. 
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